function g=prox_abs(x,a,b,c)
[num,~]=size(b);
xi=unidrnd(num,1);

a1=a(xi,:);
b1=b(xi);
z=norm(a1,2)^2;
g = x+(prx(a1*x-b1,c*z) - (a1*x-b1))*a1'/z;

function y=prx(x,c)
y=0;

if x+c<0
    y=x+c;
elseif x-c>0
    y=x-c;
else
    y=x;
end

        


